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The entanglement entropy of two gapless non-interacting fermion subsystems is computed ap- 
proximately in a way that avoids the introduction of replicas and a geometric interpretation of the 
reduced density matrix. We exploit the similarity between the Schmidt basis wavefunction and 
superfiuid BCS wavefunction and compute the entropy using the BCS approximation. Within this 
analogy, the Cooper pairs are particle-hole pairs straddling the boundary and the effective interac- 
tion between them is induced by the projection of the Hilbert space onto the incomplete Schmidt 
basis. The resulting singular interaction may be thought of as "lifting" the degeneracy of the single 
particle distribution function. For two coupled fermion systems of linear size L, we solve the BCS 
gap equation approximately to find the entropy S w (w 2 /t 2 ) logL where w is the hopping amplitude 
at the boundary of the subsystem and 2t is the bandwidth. We further interpret this result based 
upon the relationship between entanglement spectrum, entropy and number fluctuations. 

PACS numbers: 71.10.-w, 03.67.-a 



INTRODUCTION 

Quantum entanglement seems to provide an important 
connection between several distinct fields of physics rang- 
ing from conformal field theory [l[ and topologically or- 
dered phases in quantum field theories [2| to quantum 
gravity Q . One important quantity is the entanglement 
entropy, computed for a finite subregion of a quantum 
field theory (QFT) or many body system. If, in a QFT in 
d spatial dimensions, a distinguished region A of volume 
L d is formed, it follows that the degrees of freedom which 
reside exclusively in the region A will appear to be in a 
mixed state. The degree of mixing may be characterized 
by the entanglement entropy, S = — trpln p, where the re- 
duced density matrix p = tr^^ 1 0) (0 1 has been formed by 
tracing over the degrees of freedom of the ground state, 
|0), exterior to the region A. 

Entanglement entropy typically obeys an area Jaw 
and is proportional to the area of the bounding sur- 
face (L d ^ 1 ), although several variants are possible de- 
pending on the underlying particle statistics (fermion or 
boson) and dimensionality. 1 + 1-dimensional CFTs — 
which describe critical spin chains, Luttinger liquids and 
other massless theories — have pointlike bounding sur- 
faces; however, the entanglement entropy was shown to 
depend universally upon the central charge of the the- 
ory and to diverge logarithmically with the length of 
the distinguished region [3, [H| . Specifically, the entropy 
for an open chain is given by a universal expression, 
S = |lnL/e where c is the central charge and e is a 
spatial cut-off. 

Entanglement related quantities have also become a 
way of characterizing strongly correlated and topologic 
phases in condensed matter physics. Entanglement en- 
tropy diverges at a critical point and has been used to 
identify critical phases (f|. Topologic entanglement en- 



tropy (a subdominant term in the entropy) and the en- 
tanglement spectrum (the complete set of eigenvalues of 
the reduced density matrix) have been shown carry in- 
formation about the statistics of the underlying quasi- 
particles in topologic phases [7|. Recently, the entangle- 
ment spectrum of free fermion models has been linked 
to the spectrum of edge modes in topologic insulators 
Lastly, important connections between noise, full- 
counting statistics and entanglement in free fermion mod- 
els has emerged [9j. 

The universal expression for entanglement entropy of 
a CFT has a purely geometric origin. As originally dis- 
covered by Callan and Wilczek [i[ (based upon work by 
Cardy and Peschel [l(|) entanglement entropy, computed 
through a kind of replica trick, is equivalent to thermal 
entropy of a CFT mapped to a conical manifold with a 
periodic imaginary time coordinate. The appearance of 
central charge in the expression for the entropy is a conse- 
quence of the conformal anomaly and the curvature sin- 
gularity at the cone apex. The log L entropy dependence 
may be thought of as an extensive factor arising from 
the "volume" of a one dimensional gas that is quantized 
along radial equal time slices. In effect, radial degrees of 
freedom are quantized in a box of size log L. 

But there are several exceptions to the "geometry- 
entropy" connection — that is, 1-d systems exhibiting en- 
tropy proportional to log L but lacking any obvious geo- 
metrical constructions for the entropy (specifically they 
lack translation invariance or even expression as a con- 
tinuum QFT; see the discussion in reference [H, section 
3.3.) The most stri king example is found in the work 
of Refael and Moore 11| who examined entanglement in 
spin chains at the infinite disorder fixed point. In this 
limit, the wavefunction is hierarchically organized into 
products of maximally entangled RVB-likc singlets, and 
the entropy may be computed by enumerating the pairs 
that straddle the boundary between subsystems. For the 
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disordered XY model (disordered fermions), the entropy 
is found to be § In 2 In L. 

This raises the question of whether it is possible to 
generally understand the logarithmic behavior for 1-d en- 
tropy by enumeration of entangled pairs (of something). 
In this manuscript we explore building entanglements 
by weakly coupling two fermion systems, and enumer- 
ating the entangled pairs of fermions as a function of 
coupling. In particular, for two coupled gapless fermi 
systems (see figure 1), states within an energy w (the 
coupling strength) of the fermi surface should be hy- 
bridized and contribute to the entropy. This approach 
suggests expressing the entropy as a perturbation series 
in w. However, the remarkable exact solution by Eisler 
and Peschel [12j shows that this is impossible — the en- 
tropy is not analytic about w = 0, preventing a simple 
perturbative expansion. In this manuscript we suggest 
that the appropriate description at weak coupling is a 
collective, many-body effect similar to a phase transi- 
tion, and to compute the entanglement entropy we adopt 
a procedure analogous to a BCS mean field calculation. 
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FIG. 1: Two coupled fermi chains. Fermions hop between 
sites with hopping amplitude, t, and hop between chains with 
hopping amplitude, w. 

To understand the origin of the weak coupling singu- 
larity, consider two weakly coupled fermi systems as de- 
picted in figure 1. If the fermi systems are noninteract- 
ing, Peschel [1 51 ] has introduced a way of computing the 
reduced density matrix, pa, of subsystem A. Following 
15], pa has the form of a free fermion density matrix 
with a modified lattice kinetic energy 

p A = -|e- ( i) 

where c x (cj.) destroys (creates) fermions in subsytem A 
and obeys the conventional fermion algebra. The eigen- 
values of the matrix K xy may, in turn, be computed from 
the eigenvalues of the ground state free fermion correla- 
tion matrix 

n xy = (c y cl) (2) 

with x, y restricted to subsystem A. Denoting the eigen- 
values of n xy by {rife}, the eigenvalues of K xy are: 

tffc=]og(l/n fe -l) (3) 



Figure 2 is a depiction of the spectrum of for w ^ 
and w — 0. The form of rik is similar to the fermi func- 
tion at finite temperature (w ^ 0) and zero temperature 
(w = 0), although the index, k, is dimensionless. The en- 
tanglement entropy, S, is proportional to the width of the 
distribution rik at the pseudo fermi surface of one subsys- 
tem (fop = L/2 for L fermions in 2L total sites) in anal- 
ogy to thermal entropy at finite temperature. If the sub- 
sytems A and B are disconnected, n xy is formed from the 
unperturbed eigcnfunctions of the (disconnected) subsys- 
tem A: 

L/2 

nxy = J2 4>i(x)<i>i(y) (4) 

m— 1 

The eigcnfunctions of n xy are simply the complete set of 
unperturbed free particle eigenfunctions {4>m( x )} with 
eigenvalues (for m > L/2) or 1 (for m < L/2). Conse- 
quently, the eigenvalues of K xy are all singular, leading 
to the non-analytic behavior of the density matrix (fTJ) . 

The key to computing the entropy is then to describe 
how the extensive 0(L) degeneracy in the eigenvalues 
of n xy — in effect, the entanglement spectrum — is lifted 
for nonzero w. We will argue that the degeneracy is 
resolved in a manner analogous to the weak coupling 
BCS superconducting instability. Although the original 
fermions are noninteracting, the projective constraint of 
the Schmidt basis introduces an effective interaction be- 
tween particles and holes belonging to different subsys- 
tems. Since this pair interaction is singular in energy, a 
solution to the BCS gap equation results in a gap with a 
logarithmic behavior: A ~ 0(\ogL/ L). 
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FIG. 2: Depiction of eigenvalues of correlation function ma- 
trix (eq. [2]) for an L — 100 lattice. For w = 0, the distribution 
is a step function (dotted line). For w ^ 0, the width of the 
distribution (shown greatly exaggerated) is proportional to 
the entanglement entropy. 

The goal of this manuscript is to compute the logL 
entropy common to many 1-d systems without resorting 
to replicas and geometric interpretation of the partition 
function. To accomplish this we exploit the similarity 
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between the Schmidt wavefunction and superfiuid BCS 
wavefunction, where the "glue" that stabilizes the BCS 
wavefunction is a type of two particle forward scatter- 
ing across the weak link. The weak link hamiltonian 
is introduced to organize perturbation theory leading to 
an effective hamiltonian, which is then solved within a 
mean field BCS ansatz. Although our solution accom- 
plishes the stated goal — exhibiting the log L entropy — it 
fails in another way, when compared with exact and nu- 
merical results for entanglement entropy. The entropy of 
fermionic systems with a weak link has been studied in 
many ways culminating in the tour-de-force 

exact solution by Eisler and Peschel [HI for noninteract- 
ing fermion chains with arbitrary coupling strengths. For 
noninteracting chains, the exact solution shows that the 
entropy remains logarithmic in L even for weakest cou- 
plings explored. Early numerical work suggested that 
entropy was approximately quadratic in the weak link 
strength — 0(w 2 ) for w << t — but the exact solution ex- 
hibited a surprising non-analytic feature, 



(w 2 /t 2 ) log (t/w) logL 



(5) 



in this limit. In contrast, our mean field solution, which 
relies on imposing a uniform (momentum independent) 
gap, yields an entropy, S « {w 2 Jt 2 )logL. In our conclu- 
sion, we speculate on the cause of this discrepancy. 



It has been known for some time that a system of free 
fields (bosonic fields or fermionic fields) partitioned into 
two nonintersecting spatial subsystems has a coherent 
state wavefunction that resembles the BCS wavefunction 
of superconductivity. This wavefunction dates back to 
earliest explorations quantum fields in curved space and 
black hole quantum mechanics by Fulling, Parker, Un- 
ruh and Hawking [l9j . Klich and others have indepen- 
dently developed this type of wavefunction in connection 
with entanglement entropy and noise in condensed mat- 
ter systems I20l \21 \, and the density matrix renomaliza- 
tion group [22j. 

First consider the hamiltonian ([6]) with w = t, a single, 
contiguous, tight binding fermion system with length 2L. 
We will denote by <j> m (x) the eigenfunction of the correla- 
tion function, n xyi with x and y restricted to subsystem 
A and eigenvalue n m : 



n xy (j)m(y) = n m (j) m {x) 



(8) 



Klich |20( has shown that the groundstate wavefunc- 
tion may be written in the Schmidt form: 



m = n 



— f A\ f B 

L mJ m J L- 



Ma\0)b (9) 



THE SCHMIDT BASIS FOR FERMIONS 

In this manuscript, we study the entanglement entropy 
of a 1-d noninteracting fermion system (A) connected by 
a weak link to a second identical system (B). The model 
hamiltonian is: 



H = -t 



/ at q i ( * t 1 ■, 



(i,j)\a=A,B 
I At B , B\ A\ 



(6) 



where i and j are site indices and A and B denote the 
two identical systems coupled through a hopping ampli- 
tude, w. Each subsystem consists of L sites and is taken 
to have fixed boundary conditions at its ends. Further- 
more, we will restrict our calculations to the case of L 
fermions in 2L sites. This 1-d model may also be thought 
of as the effective 1-d model obtained by coupling two 
higher dimensional models by a weak link and identify- 
ing left /right moving modes with in/out s-wave scatter- 
ing states 13J, [l4[ . The Schmidt decomposition theorem 
guarantees that the groundstate of the hamiltonian (|6]) 
can be written as a sum with the form: 



a a \a) A \a) B 



(7) 



where \ct) a, for instance, are specially chosen states (the 
Schmidt basis) consisting of linear combinations of states 
drawn exclusively from subsystem A. 



where the fermion operator, fm creates a fermion 
in mode (j) m (x) in subsystem A (B); the vacuum states 
\0)a and |0)b are the zero particle state of subsystem A 
and the completely filled state of subsystem B, respec- 
tively. Thus the groundstate may be thought of as a 
BCS condensate of particle hole pairs, but with free par- 
ticle wavefunctions that are the eigenstates of n xy . From 
this wavefunction, it can shown that the entanglement 
entropy is proportional to the width of the distribution 

Based on the similarity between equation ([9]) and the 
BCS wavefunction, we attempt to compute the entangle- 
ment entropy from an appropriate BCS, mean-field start- 
ing point. For a 1-d fermion chain, this would amount 
to showing that the distribution, n m , akin to the su- 
perconductivity coherence factor, \vk\ 2 , has a width that 
depends logarithmically on the system size, L. 

There is, of course, a significant difference between the 
interacting fermions in the BCS superfiuid and the sys- 
tem of free fermions described by the hamiltonian ([6]). 
Attractive interactions between fermions in the super- 
fluid destabilize the fermi sea and lead to the BCS coher- 
ent state. Before turning to the question of the pairing 
"glue" that stabilizes this wavefunction, we mention a 
peculiarity of the Schmidt basis that is suggestive of an 
effective interaction between pairs. Since the correlation 
function, n xy , depends implicitly upon the total number 
of fermions, N, the Schmidt basis is a representation of 
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free fermions only for a specific TV. Consider the hamil- 
tonian (J6j> transformed to the Schmidt basis through: 



(10) 



x£A 



The free particle wavefunctions, 4>i(x), and eigenvalues, 
ni, also depend implicitly upon N; it follows that the ki- 
netic energy matrix of the transformed hamiltonian de- 
pends upon N. Within this representation, adding or 
removing a fermion shifts the energies of all of the other 
fermions — an effect we typically ascribe to interactions 
between fermions. 



PROJECTION ONTO SCHMIDT BASIS AND 
THE PAIRING INTERACTION 

In this manuscript, we argue that an effective particle- 
hole pairing interaction stabilizing the wavefunction @ 
arises from projective constraints imposed upon the free 
fermion hilbert space. There are many examples of ef- 
fective interactions arising this way, leading to some of 
the most interesting phenomena in condensed matter 
physics: The fractional quantum Hall effect, resonating 
valence bond antiferromagnetism and superexchange. In 
the latter case, strong coulomb repulsion removes from 
the hilbert space states in which two fermions (of oppo- 
site spin) occupy the same atomic orbital. The effective 
low energy theory is the isotropic Hcisenberg interaction 
between the spin degrees of freedom. In some sense, free 
particle kinetic energy and the projective constraint of 
no double occupancy have been " traded in" for an inter- 
action between particles. 

Our central result is that the logarithmically divergent 
entanglement entropy of the 1-d fermion system (j6|) may 
be derived from a weak coupling limit (w ~ 0) by intro- 
ducing a constraint on the hilbert space of the combined 
A and B subsystems; namely, it is forced to conform to 
the many-body Schmidt basis. The importance of the 
Schmidt basis for entanglement is that each noninteract- 
ing many body fermion state for subsystem A appear- 
ing in the wavefunction (a product of / t's) is correlated 
with exactly one unique complementary fermion state for 
subsystem B (a product of f Bjf, s). 

When the two systems are decoupled (w = 0), the 
Schmidt state is trivially: 



L/2 

n cfcf io> A |o> fl 
fc=i 



where 



L 



x=l 



L + l x 



(11) 



(12) 



(Note that the condition fg = c£ may be chosen when 
the systems are decoupled.) Equation II II corresponds to 



the first state depicted in figure [3] Since the state (fTTj) 
continues to contribute to ground state when w ^ 0, the 
next contributions must be states of the type depicted 
in the second and third states of figure [3] For instance, 
the first state depicted in figure 0] cannot contribute to 
the ground state because the zero particle-hole state for 
system B already appears in the ground state, uniquely 
correlated with the zero particle-hole state for system 
A. Similarly, the second state in figure [3] precludes the 
second state in figure 01 and so on. 
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FIG. 3: States allowed in the Schmidt basis for w ^ 0. The 
presence of the first state (N A = N B = 0) forces N A - N B to 
be even for all other states contributing to the groundstate, 
where N a is the number of particle-hole pairs in system a 
{=A,B) 



Denoting the number of particle-hole pairs in system a 
(= A, B) by N a , the projective constraint on the Hilbert 
space is 

N A -N B 



even integer 



(13) 



Of course, the total number of particles must also be 
conserved. 

This constraint is most naturally expressed kinemat- 
ically. Figure [5] depicts two possible processes allowed 
by the weak link term of the hamiltonian (see eqn. (fT5|) 
below): the scattering of a positive energy particle in A 
into a positive energy state in B, and the scattering of 
a negative energy hole in A into negative energy state 
in B. If these two processes (and their time reversed 
counterparts) are blocked, scattering occurs only within 
the constrained subspace. Next we consider how to make 
an effective hamiltonian that is consistent with the con- 
strained hilbert space. 

Following a technique introduced by Anderson for the 
renormalization group in the Kondo problem 23j , an in- 
teraction between particles is introduced that preserves 
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FIG. 4: States not allowed in Schmidt basis. The first state 
is not allowed because N B = already appears uniquely cor- 
related with N A — in figure Similarly, the second state is 
not allowed because the state depicted for B (also an N B = 
state) already appears correlated with a different state of A. 



the Dyson equation for the many-body T-matrix in the 
projected hilbert space. 
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FIG. 5: The constraint, N A - N B even, may be implemented 
kinematically by blocking the scattering of positive (negative) 
energy states in A into positive (negative) energy states in B, 
and vice versa. 



satisfies a Dyson equation: 



T(w) = H W + H w G {lo)T{lo) 



(15) 



Following |23| we try to find a new interaction operator, 
V, replacing the weak link hamiltonian, H w , so that the 
Dyson equation is satisfied in a projected hilbert space 
satisfying the constraint (|13|) . Consider the projection 
operator, P, which projects a state onto the subspace 
satisfying (|13|) . and with a fixed total number of parti- 
cles equal to L. Denote by O an operator O that has 
been projected onto this subspace; that is, O = POP. It 
is shown in [23[ that the T-matrix restricted to the pro- 
jected subspace satisfies a Dyson equation also restricted 
to the subspace, 

T(w) = V + VG a {u)T(w) (16) 

where the effective interaction, to 0(w 2 ), is: 



V = PH W P + PH W (1 - P)G a H w P - 



(17) 



The first term of equation p7|) describes the creation 
of a particle hole pair that straddles the A/B subsystems. 
The second term of (fT7|) describes fluctuations out of the 
projected subspace. Singular behavior in energy of the 
second (fluctuation) term in the effective interaction V 
is the feature that ultimately leads to the logarithmic 
dependence of the entanglement entropy on L. 
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FIG. 6: A-B particle-hole pair at [L — m, m) scatters into a 
state (m', L — m'). The intermediate state consists of a parti- 
cle hole pair, (m', L — m), in A, existing out of the projected 
subspace. 



The details of the projective RG are given in [23j; here 
we aim to give only a brief summary as it applies to 
the present case (which does not involve RG.) The T- 
matrix, T(u>), is one way of linking the bare green's func- 
tion Gq(oS) to the exact one, G(to), and is defined by the 
following equation: 



G(u) = GoM + Go(w)T( w )G (w) 



(14) 



Here, the green's function is taken to be the full, many- 
body, resolvent operator: G _1 = ui — H. The T-matrix 



To analyze this effective interaction, we rewrite H w in 
momentum space: 



H,„ = -w- 



- M mM :(c^cg,+c^cf n ,) (18) 



where 



c A ' B (x) 



sm rn c ™ (19) 
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and M„ 



sm -rrr sin -tVt ■ The kinetic energy, Hq, 



L+l 



L+l ■ 



Hq — e m ( 



c A ^c A +c B ^c B ) 



(20) 



m— 1 



-2* COSfiZL. 



where e m 

We consider a situation, as depicted in figure ©, where 
an particle-hole pair at (L — to, to) scatters into 

a state (m',L — to'). The intermediate state consists 
of a particle hole pair, (to',L — to), in A, existing out 
of the projected subspace. Restricting the scattering as 
described, the relevant term in the effective interaction 
(PU is: 



V 



Aw 2 



r M C B r M r B 



m,m' >0 



w - (e m + Em') 



V V ,c A 



A\ B 
L—rn' ^m' 



(21) 



This interaction, as in the phonon mediated BCS-theory, 
is a frequency dependent retarded interaction, and we 
concentrate on the low frequency behavior uj s» where 
the interaction is attractive. 



BCS SOLUTION OF THE EFFECTIVE 
HAMILTONIAN 

When the two subsystems are disconnected (w = 0), 
the wavefunction © degenerates to two disconnected 
fermi seas as expected. The eigenfunctions of n xy , as 
described above in ©, become the free particle eigen- 
functions of the half-chain, and n m = 9(L/2 — to). The 
wavefunction (|9]) then becomes 



l^> = II(VT" 



cf-m)|0)A|0)i 



(22) 



the variational energy becomes: 

U = E ( 2e ™"'» ~ A m -y/n m (l - n m )) (24) 



Minimizing the energy, U , with respect to n m and solving 
for n m gives coherence factors that take the usual BCS 
form: 



9 TP ' 



(25) 



where E 2 n = A 2 n + (2e m ) 2 . 

Next we turn to solving the gap equation (f23|). Al- 
though the scattering amplitude, V m . m i, is an energy de- 
pendent interaction, we will seek an approximate uniform 
solution to the gap equation. V m ,m' has a maximum value 
when to, to' ss L/2; that is, when scattering events are 
close to the fermi points of subsystems A and B. We eval- 
uate the scattering amplitude in the static limit, uj ~ 0, 
and for a low energy initial pair state (e m ~ 0). Lineariz- 
ing about the fermi point, to = L/2, the kinetic energy 
is 



2irt 



-p p = ±l,±2... 



and the effective interaction, V m m ', becomes: 



V m 



W 2 1 

Lt to' 



(26) 



(27) 



where L+l has been replaced by L and a numerical 
factor has been absorbed into w. 

Following the standard BCS procedure and choosing a 
uniform solution for the gap, A = A m , the gap equation 
(l2lif becomes: 



1 w 2 1 1 
2Ti ^ mil, 



(28) 



where c n 



4 



create and destroy simple band fermions 



in the disconnected subsystems. For weak coupling, 
w ~ 0, the pairing approximation of BCS theory as- 
sumes that only the coherence factors y/T— n m and y/n m 
in (|22l) — these are the familiar Uk and Vk coefficients of 
superconductivity — are modified from their discontinu- 
ous distributions at w — 0, and the wavefunctions corre- 
sponding to c m , c m are not. 

Following the BCS procedure, we take the wavefunc- 
tion (|22l) as a variational wavefunction for the hamilto- 
nian, H = Hq + V . To simplify the notation, constants 
and the singular part of the effective interaction have 
been grouped into the scattering amplitude V m , m > on the 
second line of (|21[) . Defining a "gap" function, 



A m = 



Vrn.rn' \f^rn' (1 ^hn' ) ) 



(23) 



To solve the gap equation, we convert the sum (|28| to an 
integral and cut off the energy integration at a scale, at, 
where a(< 1) is a numerical constant, giving: 



i - i^L 

~ 4 Lt 



dc 



t /L e ^e 2 + (A/2) 2 



(29) 



Looking for a self-consistent solution, we approximate the 
integrand by 2/(Ae) and integrate to get: 



A 



2,11 



In L + c(a) 



(30) 



where c(a) is a constant depending upon the cut-off. This 
solution is self-consistent in that A is comparable to the 
energy cut-off (justifying truncation of the integrand) for 
large L. 
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Within the Sommerfeld approximation, the entropy is 
proportional to the product of the gap and the density of 
states at the fermi point, N = L/2irt, and we arrive at: 



S « ANn 



w 



InL 



(31) 



This result agrees qualitatively with numerical computa- 
tions of the entanglement entropy for 1-d fermions in a 
weak link geometry (l6j . Specifically, the entropy is log- 
arithmic and the prefactor is approximately quadratic in 
w. Other numerical computations of entropy for 2-d and 
3-d systems with weak links (i.e. quasi-one dimensional 
systems) have also exhibited approximately quadratic be- 
havior [13|. However, the exact solution (and earlier 
work) informs us that the prefactor is non-analytic and 
approximately proportional to w 2 In (1/w) [Hf. Even in 
the original numerical work fl6j j . Peschel noted a very 
small, persistent deviation from quadratic dependence 
upon the prefactor. We return to this point in the Dis- 
cussion section. 



NUMBER FLUCTUATIONS AND ENTROPY 

In a 1-d translation invariant fermion system, it has 
been rigorously established that the entanglement en- 
tropy of a subsystem is proportional to the number fluc- 
tuations in that subsystem 0, [2^, 21, 24 1. For such a 
continuum 1-d system, the entanglement entropy of a 
length L domain is proportional to logX with a univer- 
sal prefactor that depends upon the central charge of the 
underlying CFT. For subsystems separated by a bond or 
site defect, the entropy remains logarithmic but with a 
modified pre-factor [12j . 

However, in such systems lacking translation in- 
variance the relationship between entropy and num- 
ber fluctuations — more generally, the complete set of 
cumulants — becomes complicated. Klich and Levitov 
studied entanglement entropy in a gated quantum point 
contact (QPC) and have shown that entropy may be ex- 
pressed as an asymptotic series in the cumulants (see also 
reference [25[ for refinements.) For a QPC with perfect 
transmission, all but the first and second cumulants van- 
ish, and entropy is simply proportional to number fluctu- 
ations. For barriers with nonzero reflection coefficients, 
relevant to the weak link geometry, the FCS generating 
function is non-Gaussian and all even cumulants gener- 
ally contribute to the entropy. 

In this section we explore the number fluctuation- 
entropy connection based upon perturbation theory ap- 
plied to weakly coupled fermion chains. The number fluc- 
tuations in one subsystem SN 2 = (N 2 ) — (N) 2 is propor- 
tional to log L, but with a prefactor that is analytic as 
w — > (a result derived earlier in [26jj). Since the en- 
tropy in this limit is known through the exact solution 
to depend nonanalytically on w, number fluctuations and 



entropy cannot be proportional. The ratio of entropy to 
number fluctuations turns out to be 



S , t 
— — t oc In — 

SN 2 w 



(32) 



and is divergent as w — > 0. This result suggests that the 
nonanalytic behavior of the entropy can only be captured 
by resumming the logarithmic terms in the cumulant se- 
ries for the entropy to all orders in w. 

Fluctuations in N may be computed perturbatively in 
w using straightforward canonical methods. Defining the 
density-density correlation function 



n(x,y,t) = (Tp(x,t)p(y,0)) 



(33) 



where p(x, t) = c^(x, t)c A (x, t), the number fluctuations 
may be computed from 



L 

(N 2 ) = n(z,y,o) 

x,y=l 



(34) 



Denoting by G the single particle propagator for system 
A (not to be confused with the many-body resolvent op- 
erator G(w)), G may be written in the interaction pic- 
ture, 

G(x,y,t) = -i(Tc A {x 1 t)c A \y 1 0)e- l S-™ H ™ (t ' )dt ') (35) 

where T is the time-ordering symbol and (. . .) refers to 
average taken in the ground states of the disconnected 
subsystems A and B. Expanding the propagator to 
0(w 2 ), 

w 2 f 00 f°° 
G 2 (x,y,0-) = -— / dh / dt 2 G (xA,-t 1 ) 

* J-oo J—oo 

x G (l,2/,t 2 )Go(l,l,ii-t 2 ) (36) 

where Go is the unperturbed propagator 

G (x, y, t) = -i{Tc A (x, t)c A \y, 0)) (37) 

Inserting the expansion of G up to second order into n, 
we obtain 

U{x,y) = Ga(x,y)G (x,y) + Go(x,y)G 2 (x,y) + ... (38) 

where we have dropped the time argument, t = 0. The 
covariance to 0(w 2 ) is then: 

sn 2 = (n 2 ) - (ao 2 - G °( x ' y) G2 ( x > y) ( 39 ) 

Expressing G 2 in spectral form, the covariance may be 
written: 



L/2 L 



5N 2 = 



M 2 , 

m.rn 



=w=i ^ m £m ) 



(40) 
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where S m = sgn(e m ) and the resonant terms m = rn! are 
excluded. The sum (|40|) may be turned into an integral 
and evaluated to give 

2 

W 

6N 2 = C^ r kiL (41) 

Note that this result is qualitatively the same as the log- 
arithmic behavior found for the number fluctuations in 
an L site subsystem of a translationally invariant chain. 
Using (|4"Tj) and ([5]) to compute the ratio of entropy to 
number fluctuation establishes the result (|32p. 

We note that fluctuations in the number of fermions, 
N, in subsystem A may also be computed from the 
Schmidt wavefunction Q within our projective ap- 
proach. Transforming 

<^ 2 > = E WcfcX^W (42) 

x.yeA 

to the Schmidt basis using the canonical transformation 
(JTU]), the fermion number covariance is found to be: 

L 

(N 2 ) - (TV) 2 - £ n m (l - n m ) « AN (43) 

m— 1 

Again, within the Sommerfeld approximation, the co- 
variance is proportional to the width of the distribution 
n m (l — n m ) appearing in the sum above, and gives the 
same result as found for the entanglement entropy in 
(I3ip . Thus within our approach, entropy and number 
fluctuations are simply proportional. 

Lastly, we return to the perturbative number fluctua- 
tion results, equations (PffJl) and (|4"Tj) and interpret them 
in a different way. The projective approach described 
previously was motivated by the observation that the ex- 
tensive degeneracy in the eigenvalues of the correlation 
function (they are all or 1) is lifted for any non-zero w. 
From the shift in degenerate eigenvalues of the correlation 
matrix (from or 1), the entropy might be computed. 
Since G2 is the first nonvanishing correction to the cor- 
relation function, the correction to the mth eigenvalue of 
the correlation function, 5n m , is found from degenerate 
perturbation theory 

6n m ^J2^)G2(x,y)^i(x) (44) 

xy 

Note that Yl,m=i & n m is the same as the expression for the 
number fluctuations, (f3U)) , and leads to the same results, 
equations (|40j) and (|4T|) . 

Thus the perturbative expression for number fluctua- 
tions is really the same calculation as degenerate pertur- 
bation theory applied to the extensive degeneracy of the 
correlation function eigenvalues for disconnected chains. 
However, inspection of the numerical eigenvalues, n m , 
shows that the perturbative corrections to n m are not 



even qualitatively correct, even though the final expres- 
sion for entropy is correct. According to perturbation 
theory, 6n rn cx 1/m for large m, whereas the behavior 
for 5n m seen numerically — and by the BCS approach 
described in this manuscript — is a "gapped" behavior, 
where n m is given by the expression (|2"5)) . 

If higher order cumulants (N ) c , (N 6 ) c etc. are also 
analytic functions of w about w = 0, it seems that the 
distribution n m can only be obtained by resumming the 
cumulant scries to all orders in w. Presumably each rel- 
evant term is logarithmically divergent in L. 

DISCUSSION 

The motivation of this work was finding an alterna- 
tive explanation for the log L entanglement entropy — one 
that does not rely upon replicas and a geometric inter- 
pretation of the resulting partition function. In the latter 
approach, entanglement entropy appears as the fictitious 
thermal entropy associated with the partition function 
evaluated in a conical geometry (or, equivalently in terms 
of corner transfer matrices.) In a conical geometry the 
log L behavior may be traced to the log-periodic coordi- 
nates that diagonalize the free particle kinetic energy in 
polar coordinates. Schematically, the entropy is evalu- 
ated for a partition function, 

Z = J D(f>e' Sa (45) 

that follows from the free particle action on a cone with 
completion angle f3 (the inverse temperature): 

S a = [ d6 [ dr((d r cf>) 2 + \(d g 0) 2 ) (46) 
Jo Jo r 

The action is diagonalized by 

</>{r,6) ~ e im< V fe ^ (47) 

where L is the radial length and e is a radial cut-off. 

The divergent log L entropy appears as a consequence 
of each momentum eigenvalue collapsing inversely pro- 
portional to its quantization "length," logL/e. Just 
as entropy in a massless quantum fluid at temperature 
is S w L/(3 — the aspect ratio of length and imag- 
inary time dimensions — entanglement entropy becomes 
S^^logL/e. 

In our calculation, entanglement entropy is computed 
from the number of particle hole pairs that straddle the 
boundary between subsystems. The particle hole pairs 
have an effective attraction and form a BCS type con- 
densate; the product of the energy gap, A, and the den- 
sity of states at the fermi level, A^o, is a measure of the 
number of participating pairs and, therefore, the entan- 
glement entropy. Within our approximations, the gap, 
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A fa jj InL, and the density of states No « L/t, giving 
our central result: 

2 

S«^-lnL (48) 

Although our calculation produces the logarithmic en- 
tropy in a novel way, it only agrees approximately with 
the exact results of entanglement entropy in a fermion 
chain. Specifically, our result fails to produce the cor- 
rect prefactor in the entropy, the exact result being 
S ss (w 2 /t 2 ) log (t/w) logL in the limit w « t. 

To understand this discrepancy we first return to the 
exact solution at the translationally invariant point (w — 
t) and consider the reduced density matrix ([1} and ef- 
fective single particle eigenvalues, K k (equation ^) de- 
scribing subsystem A. In the translation invariant limit 
w = t, the single particle eigenvalues Kk appearing in the 
quasi thermal fermi distribution (adopting the notation 
of reference mav be thought of as the quotient of 
an effective linearized energy, Uk ntk/L and an effec- 
tive temperature T — irtL^ 1 InL. (That is, Kk = 0Jk/T.) 
These results follow from the asymptotics of Peschel [28[ 
and are convincingly demonstrated in the numerical work 
of Eisler et al (see figure 1 of ref. [HI). The entropy is 
thus the product of the temperature, T, and density of 
states, L/irb, giving S ss logL. Our result, although 
arrived at in a completely different way (and in a com- 
pletely different limit, w << t), resembles this result by 
identifying the effective temperature, T, with our gap, A. 
We also note that BCS type distribution, rife, with gap 
A in our result, and a fermi distribution with effective 
temperature T = A/2 agree to 0(ef./A 2 ) and thus give 
the same result for the entropy in the limit L — > oo. The 
two distributions yield identical entropies in that that 
their widths, which determine the entropy, both vanish 
as 0(L _1 logL). In summary, our result for the distribu- 
tion rife, derived in weak coupling w « t, resembles the 
rife obtained in the exact solution for w = t. 

The discrepancy of our result with the exact result at 
weak coupling (w << t) may be traced to a qualitative 
change in the effective single particle eigenvalues, Kk, 
when w << t. In the exact solution, Wfe develops a gap 

approximately of the form u>k = \J A 2 + e|, for w « 

t. The gap A (distinct from our coherence gap, A), is 
approximately A ss hi (t/w). This is the origin of the 
In (t/w) appearing in the prefactor of the entropy. In the 
BCS approach, the distribution rife (comparable to the 
fermi distribution in the latter analysis), involves only 
the free particle kinetic energies of the lattice fermions 
which are, of course, ungapped. Although this appears to 
be a fundamental limitation of our calculation, perhaps 
an improved variational approach might capture the gap 
feature. 

To summarize these relationships: the computation of 
both number fluctuations and entropy based upon the 



Schmidt wavefunction amount to evaluating the width 
of the distribution rife. In the exact solution, when 
w = t, rife has the form of a Fermi distribution with a 
linear single particle dispersion resembling that of free 
fermions and an effective temperature T — 7rfL _1 lnL. 
Within the projective/BCS type approach described in 
this manuscript, rife has a gapped behavior given by 
equation (|25|) with an approximate linear single parti- 
cle dispersion given by the original lattice fermions and 
a gap, A sa (w 2 /t)L~ x InL. Thus these two results qual- 
itatively agree as w — > i, given the vanishing width of 
rife in both cases as L — > oo. On the other hand, the 
computation of number fluctuations (or entropy) from 
canonical perturbation theory leads to a result propor- 
tional to logL, although nk has a qualitatively different 
behavior — algebraic rather than gapped. Our projective 
approach is then an improvement over perturbation the- 
ory, and in qualitative agreement with rife in the trans- 
lation invariant case, but, as discussed in detail above, 
fails to agree with the exact solution when w << t. 

One might argue that BCS theory of superfluidity is 
mean field theory and is only expected to agree approx- 
imately, especially in the case of a one dimensional ex- 
act result. Specifically, BCS captures a thermodynamic 
phase transition, where the order parameter (the phase of 
the wavefunction representing the condensate) becomes 
fixed, stabilized by the free fluctuations in the conjugate 
variable, TV, (the number of particles in a patch of size 
equal to the coherence length.) The BCS approximation 
fails badly in one dimension and marginally in two di- 
mensions because phase fluctuations are too strong. Sim- 
ilarly, in a finite size sample, phase fluctuations are not 
frozen and, in fact, correspond to the zero point motion 
of a massive rotator. Since our computation is aimed 
at finding the finite size scaling of entropy, one might 
question using the mean field result where fluctuations 
of the order parameter have been frozen out. Moreover, 
our BCS approach entails a solution to the energy de- 
pendent gap equation that imposes a uniform gap; this 
also may be responsible for its failure to reproduce the 
correct nonanalytic singularity. 

The main ingredient of our result is the projective 
constraint on the underlying hilbert space — an approach 
that might be used in any system where a variational 
form of the Schmidt basis is postulated. Another advan- 
tage to this approach is the focus upon computing rife , in 
particular finding a reasonable (and physical) approxima- 
tion for this highly singular distribution. For noninteract- 
ing fermions, {rife} encodes all eigenvalues of the reduced 
density matrix and therefore the entire entanglement 
spectrum (a comparable argument exists for bosons). 
Understanding how the degeneracy in rife is generically 
resolved at weak coupling might provide a path to com- 
puting entropy in higher dimensions where specialized 
I-d methods cannot be used. It is also interesting to con- 
sider if the non-analytic feature in the entropy at weak 
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coupling survives in higher dimensions. If fermions in- 
teract in such a way that the system is still described as 
a fermi liquid, this technique might also provide an ap- 
proach for computing the entanglement spectrum in such 
systems. The approximation in this case might resemble 
the Bogolubov-de Gennes [lo[ formalism where interac- 
tions are decoupled into off-diagonal Hartree terms to 
describe superfluidity and diagonal Hartree terms to de- 
scribe residual interactions. 

Lastly, our perturbative calculation of the second cu- 
mulant suggests an investigation of the cumulant series 
for entanglement entropy. Since the second cumulant is 
analytic in w, it seems necessary to re-sum the series to 
all orders to capture the nonanalytic behavior exhibited 
by the exact solution. 

We thank the referees for pointing out the origin of 
the discrepancy with the exact solution and references 
2tj 29[ . This research was supported by grants from the 
Department of Energy DE-FG02-08ER64623— Hofstra 
University Center for Condensed Matter, Research Cor- 
poration CC6535 (GL) and NSF Grant no. 0705048 
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